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Abstract 

We summarise the recent progress of a project to deter- 
mine physical conditions in the jets of low-luminosity 
(FR I) radio galaxies. We model the jets as intrinsi- 
cally symmetrical, relativistic, decelerating flows. By 
fitting to deep, high-resolution VLA images in total in- 
tensity and linear polarization, we have been able to 
derive for the first time the three-dimensional distribu- 
tions of velocity, rest-frame emissivity and magnetic- 
field structure. We describe our results for four sources 
and outline the implications for theories of jet deceler- 
ation. We also summarise a conservation-law analysis 
of 3C 3 1 , which gives profiles of internal pressure, den- 
sity, Mach number and entrainment rate along the jets, 
again for the first time. Finally, we discuss the applica- 
bility of adiabatic models to FR I jets and identify the 
regions where distributed particle acceleration is nec- 
essary. 

1 Introduction 



It was first recognised by Fanaroff & Riley ( 1974) that 
the ways in which extragalactic jets dissipate energy 
to produce observable radiation differ for high- and 
low-luminosity sources. Jets in weak FRI sources 
are bright close to the nucleus of the parent galaxy, 
whereas those in powerful FR II sources are relatively 
faint until their terminal h ot-spots Tthere is also a de- 



pendence on environment (Led low & Owerl 199flh l. It 



rapidly became accepted that FR I jets m ust decelerate 
by entrainment o f the surro u nding IGM JfjaanL 1980; 
i BegelmanL Il982l: iBicknellL 1 1984 Il984 be Young . 
199(4 120041 : IRosen et all 1 19991: iRosen & Hardee . 
2000) or by injection of mass lost by stars within 



the jet volume JphinnevL 19831: Komissarov . 1994; 
Bowman. Leahy & Komissarovl. 1996). 



More recently, evidence has accumulated that FRI 
jets are initially relativistic and decelerate on kpc 
scales. FRI sources are thought to be the side- 
on counterparts of BL Lac objects, in which rela- 
tivistic motion on pars ec scales is well-established 
JUrrv&PadovaniL ri995\ Superluminal motions have 
been seen on mill i arcse cond scales in several FRI jets 



( Giovannini et all 2001) and on ar csecond scales in 
M87 (iBiretta. Zhou & OwenLll995h . In FRI sources, 
the lobe containing the main (brighter) jet is less de- 



polarized than the counter-jet lobe (Mor ganti et al 



1997). This can be explained if the main jet points to- 



ward the observer, suggesting that the brightnes s asyrn 



198 



metry is caused by Doppler beaming dLainf 
The a symmetry decrease s with distance from the nu- 
cleus I Laing et al. , 1999h . consistent with the idea that 
the jets decelerate. 

It has proved frustratingly difficult to quantify the basic 
properties - composition, velocity, density, pressure, 
magnetic field - of radio jets. The present paper sum- 
marises recent progress on a project which has, for the 
first time, allowed us to derive the three-dimensional 
distributions of physical parameters in FRI jets on kpc 
scales. There are three main elements: 

1. We model FRI jets as intrinsically symmetrical, 
axisymmetric, relativistic flows. We parameterise 
the geometry and the three-dimensional variations 
of velocity, emissivity and magnetic-field struc- 
ture and compute the brightness distributions in 
total intensity and linear polarization. By fitting to 
radio images, we can optimise the model parame- 
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ters. We refer to these fits as free models. They are 
empirical, and make as few assumptions as possi- 
ble about the (poorly-known) internal physics of 
the jets. Deep, high-resolution radio images in 
Stokes /, Q and U are essential. 

2. We then apply conservation of mass, momentum 
and energy to infer the variations of pressure, den- 
sity, entrainment rate and Mach number along the 
jets. This requires measurements of the external 
density and pressure profiles from X-ray observa- 
tions. 

3. Finally, we attempt to model the acceleration and 
energy-loss processes, adding images at mm, IR, 
optical and X-ray wavelengths to our radio data 
and models. So far, we have concentrated on the 
simplest process - adiabatic energy loss - in an at- 
tempt to quantify the amount and location of dis- 
tributed particle acceleration. 

Two X-ray - radio connections are vital to our argu- 
ment. Firstly, the bases of FRI jets emit X-ray syn- 
chrotron radiation, indicating continuing local particle 
acceleration. Secondly, the jets interact with and are 
confined by the hot intergalactic medium, whose den- 
sity and temperature can be estimated from its thermal 
bremsstrahlung emission. The radio data alone provide 
constraints on the kinematics of the jets. Combining 
these kinematics with the densities and pressure gradi- 
ents obtained from the X-ray data via a conservation- 
law analysis allows us to address the jet dynamics and 
to estimate Mach numbers, energy and mass fluxes and 
entrainment rates in the jets. 

We began this project by making a statistical study o f 
jet asymmetries in the B2 sample dLaing et al. , 19991) . 
The detailed models were fi rst applied to the well- 
known FRI source 3C 31 by Laing & Bridlel J2002ah 
and a conserv ation-law analys i s of th e same source 
was made by lLaing & Bridlel (20023), using Chan- 
dra observations by H ardcastle et al.ld2002h . Adiabatic 
mo dels were developed a nd applied, again to 3C 3 1 , 
by Laing & Bridlel J2004T) . More recently, the tech- 
niques have been extended to other sources and free 
mod els of 0326+39 an d 1553+24 were published by 



2 Free models 

2.1 Model procedure 

Our fundamental assumption is that the jets are intrin- 
sically symmetrical, axisymmetric, relativistic, station- 
ary flows. [The assumption of symmetry is critical to 
our approach and we consider it further in Section l2~2ll . 
The magnetic fields are assumed to be disordered on 
small scales, but anisotropic. There are both observa- 
tional and theoretical reasons to suppose that the lon- 
gitudinal field compon ent, at least, cannot be vector- 
ordered on large scales dBegelman. B landford & Rees, 
1984l : lLaingLll98ll:lLaing & Bridld.l2002ah . We return 
to the question of the nature of the toroidal field com- 
ponent in Section l3~4l 

The main elements of the modelling procedure are as 
follows: 

1. We parameterise the geometry, velocity, emis- 
sivity and field struc ture using simp l e, ana- 
lytical express i ons Eaing & Bridlel Ep02a; 
Canvin & LainslEooi) . 



2. For given model parameters, we derive model dis- 
tributions of Stokes IQU by integration along the 
line of sight, taking account of the anisotropy 
of synchrotron radiation in the rest frame of the 
emitt ing material , aberration and Doppler beam- 
ing (O3|2002|). 



Canvin & Laing (2004). In the present paper, we sum- 



marise the published results and also present models of 
the jets in NGC315. Modelling of one further source, 
3C 296 dHardcastle et all 1 19971) . is in progress. 



3. The model is then compared with deep VLA im- 
ages, using x 2 (summed over Stokes parameters) 
as a measure of goodness of fit and the parameters 
ar e optimised using the d ownhill simplex method 
of iNelder & Meadl dl965l) . We take a measure of 
the deviation of the brightness distributions from 
axisymmetry as an estimate of the "noise" (rather 
than the much smaller off-source rms) and obtain 
reduced \ 2 values in the range 1.1 - 1.7 over 1200 
- 1400 independent points, each with 3 Stokes pa- 
rameters. 

2.2 The assumption of intrinsic symmetry 

It is, of course, an assumption that the jets in most FRI 
galaxies are intrinsically symmetrical on small scales - 
and one that is difficult to prove in any individual case 
[there are also some clear exceptions, such as 0755+37 
i Bondietall EoOOll. It is also obvious that many jets 
are very far from symmetry on large scales. We con- 
tend, however, that the apparent asymmetries caused 
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by aberration and Doppler beaming dominate over in- 
trinsic and environmental differences close to the nu- 
clei in most sources. Evidence in favour of this hypoth- 
esis includes: 



1. 



The intensity ratio between main and counter jets 
initially decr e ases w ith distance from the nucleus 
( Laing et"all 1 1999b . reaching a value close to 
unity on scales ~10kpc in typical cases. Deceler- 
ation provides a natural mechanism to reduce the 
asymmetry if the flow is relativistic, whereas in- 
trinsic or environmental differences between the 
jets would be likely to produce constant or in- 
creasing brightness differences. The fact that the 
jets are very symmetrical on intermediate scales 
indicates that intrinsic differences are small. 

Jet/counter-jet intensity ratios are well corre- 
lated with other indicators of relativistic motion, 



such as the core/extended flux ratio ( Laing et al 
1999h . 



In those cases where superluminal motion is 
observed, it is always on the side of the 
brigh t er kpc-scale jet fe iretta. Zhou & QwenL 
199.4 iGiovannini et al.Ll200lh . 



The observation that the lobe containing the 
brighter kpc-scale jet is less dep olarized than the 



coun ter-jet lobe in FRI sources ( Morgandetal 
1997) indicates that the brighter jet is seen 



through less magnetoionic material and is there- 
fore on the near side, as expected. 

5. There is a characteristic asymmetry in the degree 
and direction of polarization (corrected for Fara- 
day rotation), which is correlated with that in total 
intensity. We show below that relativistic aberra- 
tion provides a natural explanation for this effect. 

Fig.^illustrates the characteristic asymmetries in total 
intensity and linear polarization for the four sources in 
our study. They are arranged in order of increasing an- 
gle to the line of sight, 9 (as deduced from the model 
fits). The asymmetries in total intensity decrease with 
distance from the nucleus and are, of course, greater 
at smaller 9. In addition, there is an asymmetry in 
the polarization structure: the main jet base shows a 
polarization minimum on-axis which extends consid- 
erably further from the nucleus than any equivalent in 
the counter-jet. Our (intrinsically symmetrical) models 



successfully reproduce this difference (see Section l2l4l 
and we believe that it, like the intensity difference, is a 
signature of relativistic aberration. 

2.3 Breaking the degeneracy between velocity and 
angle 

The key to our modelling approach is its use of asym- 
metries in linear polarization as well as total intensity 
in order to break the degeneracy between velocity and 
angle. The standard formula for the ratio of jet and 
counter-jet flux densities (assuming isotropic emission 
in the rest frame) is: 

f l+(3cos9 \ 2+a 
j/cj " \1-PcobO) 

where a is the spectral index and = v/c. There is 
therefore a degeneracy between velocity and angle to 
the line of sight. The relation between the angles to the 
line of sight in the rest frame of the flow, 9' and in the 
observed frame, 9, is: 

sin(9j = [T(l -p cos 9)]' 1 sinfl (main jet) 
sin 6' ci = [T(l + 13 cos 9)Y l sin 9 (counter-jet) 

where V = (1 — /3 2 ) _1//2 is the Lorentz factor. The 
observed polarization is in general a function of 9' if 
the field is anisotropic. If we know the field structure 
a priori, then we can solve explicitly for (3 and 9. We 
take the example of a field which is disordered on small 
scales but confined to a plane perpendicular to the jet, 
with equal rms along any direction in the plane. In this 
case, there is no variation of the degree or direction of 
polarization across the jet. For a = 1, the total and 
polarized flux densities per unit length in the emitted 
frame are: 



P' 



iv~(l + cos 2 9') 



p Q K sin 2 9' 



where K is a constant dLaingl [198(1 fl98ll The ratios 
of observed total and polarized intensity for the jet and 
counter-jet are: 

3 



~ 2 sin 2 I 



2- [T(l -/3cos0)] 
2- [r(l + /3cos#)]- 2 sin 2 i 



1 + (3 cos 9 
1 -/3cos6> 
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NGC 315 




Figure 1 : The observed total-intensity and polarization structures of the four modelled sources. The sources are 
arranged from top to bottom in order of increasing angle to the line of sight, 9, as deduced from the model fits. 
The left-hand panels show vectors with lengths proportional to the degree of polarization, p and directions along 
the apparent magnetic field superimposed on false-colour plots of total intensity. The right-hand panels show 
false-colour plots of p, in the range (blue) - 0.7 (red). Note: (a) the asymmetry in total intensity decreases with 
increasing 6; (b) there is a characteristic difference between the polarization structures of the main and counter-jets 
(see te xt) and (c) 0326+39 has no parallel-field edge, u nlike the other three sources. Data from Canvin & Laingl 
< 2004h (0326+39 and 1 553+241 Fling & Bridle! J2002ah f3C 31 ) and Cotton et al., in preparation (NGC315). 
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These equations can be solved numerically for (3 and 
6 (cf. iBondi et alJEoOflh . In general, we must fit the 
field configuration and therefore need to introduce ad- 
ditional parameters to describe its anisotropy. These 
also determine the variation of polarization transverse 
to the jet axis, and can be estimated independently of 
the velocity and angle if the jets are well resolved in 
this direction. 

The selection of sources to be modelled is determined 
by the requirement for high signal/noise ratio and good 
transverse resolution in /, Q and U in both jets. It is 
also essential to correct for Faraday rotation. 

2.4 Comparison of models and data 



The four sources we 



have modelled 
0326+39, 



are 3C31 
1553+24 



Laine & Bridle, 2002a), 

Canvin & Lainl 12004b and NGC315 (Cotton et 
al., in preparation). In Figs 12-0] we compare model 
and observed images at the same resolution. Figs|2]and 
[3] show the full modelled areas at moderate resolution 
and display both total intensity and linear polarization. 
They demonstrate that we can successfully fit the 
observations in all cases, limited primarily by small- 
scale (and often non-axisymmetric) fluctuations. The 
inferred angles to the line of sight are shown in the 
figures. Fig. 0] shows the inner parts of the model and 
observed / images at higher resolution. Here, the 
models fit the average intensities very well, but cannot 
reproduce the complex small-scale structure. 

3 Physical parameters deduced from the free 
models 

3.1 Geometry 

Our qualitative conclusions about the geometry of FR I 
jets are not new, but our estimation of angle to the line 
of sight allows us to deproject their shapes and to de- 
rive intrinsic parameters. FRI jets are initially narrow, 
expand abruptly and then recollimate. We refer to this 
part of the flow as the flaring region. The shape of 
its outer isophote can be fitted accurately by a cubic 
polynomial function of distance along the axis. Fur- 
ther from the nucleus, the outer region is conical. In all 
sources except 3C31 the opening angle is very small, 
so that the outer regions are almost cylindrical (note 
that we model only the flaring region in NGC315). 
The deprojected shapes of the sources are shown in 
Fig. |5] The very small angle to the line of sight inferred 




NGC 31 5 data 




Figure 2: Comparison between the model and ob- 
served brightness distributions for 1553+24 (upper two 
panels) and NGC 315 (bottom two panels). Vectors 
with lengths proportional to the degree of polarization 
and directions along the apparent magnetic field are su- 
perimposed on false-colour plots of total intensity. The 
resolutions are 0.75 and 2.35 arcsec FW HM, respec- 
tively. Data and models from Canvin & Laing l (2004) 
(1553+24) and Cotton et al, in preparation (NGC 315). 
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■ 553-21 data 



■ - 3 dag 



• 553-21 modal 








ft =64 deg 



Figure 3: Comparison between the model and ob- 
served brightness distributions for 3C 3 1 (upper two 
panels) and 0326+39 (bottom two panels). Vectors 
with lengths proportional to the degree of polarization 
and directions along the apparent magnetic field are su- 
perimposed on false-colour plots of total intensity. The 
resolutions are 0.75 and 0.5 arcsec FWHM. respec- 
tively. Data and models fro m Lain g & Bridle! ( 2002ab 
(3C31) and lCanvin & Laingl J2004 0)326+391. 



032B *-3B model 



Figure 4: Comparison between model and observed 
total-intensity images at high resolution. Note the non- 
axisymmetric and knotty structure observed where the 
jets first brighten significantly. This cannot be de- 
scribed in detail by our axisymmetric models, but the 
average structure is reproduced correctly. NGC315 
and 3C31 clearly show highly-collimated inner jets: 
these are probably present in the other two sources, but 
are too faint or small to image. References as in Figs|2] 
and|3l 
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for 1553+24 means that its intrinsic shape is very nar- 
row (see Section l3~6T l. 

3.2 Velocity 

Fig. |5] also shows the inferred velocity distributions. 
At the closest distance from the nucleus where we can 
model the jets, all of their on-axis velocities are close 
to f3 ~ 0.8. They all decelerate rapidly over a short 
distance in the flaring region, but at different distances 
from the nucleus. NGC315, 0326+39 and 1553+24 
then have roughly constant velocities in the range j3 m 
0.1 - 0.4 until the end of the modelled regions. 3C 31 
continues to decelerate, but more slowly. 

There is clear evidence for transverse velocity gradi- 
ents in 3C 31 and NGC 315. In the other two sources, 
the gradients are poorly determined, but in all four 
cases, the data are everywhere consistent with a con- 
stant ratio of edge/on-axis velocity sa 0.7. There are 
no obvious low-velocity wings. 

3.3 Emissivity 

The longitudinal emissivity profiles generally flatten 
with increasing distance from the nucleus. All of the 
sources show an anomalously bright region, with com- 
plex, knotty sub-structure, within the flaring region 
(e.g. Fig. |6]l. This has an abrupt start, correspond- 
ing to the flaring point of the jet JParma et al. , Il987: 



Laing et all |f999) and often also a clear end: it is as 
if an extra component of emissi vity is added to an un- 
derlying, smooth distribution. Laing & BridQ d2002ah 
suggested that the sudden brightening coincides with a 
change in jet collimation properties in 3C31, but the 
transverse resolution is marginal for that source, and 
there appears to be no obvious change in collimation 
in NGC 315, the only well-resolved case. 

The jets are all intrinsically centre-brightened, with a 
typical edge/centre emissivity ratio « 0.2. 

3.4 Field structure 

We have, for the first time, estimated the three- 
dimensional field structures of the jets rather than their 
projections on the plane of the sky. The dominant 
field component at large distances from the nucleus is 
toroidal. The longitudinal component can be signif- 
icant close to the nucleus, but decreases further out, 
sometimes becoming negligible. The behaviour of the 
radial component is peculiar, and varies from source 
to source. This field evolution is qualitatively consis- 




Figure 5: The velocity structures deduced for the four 
modelled sources. The labelled wedge indicates val- 
ues of j3 in the range - 0.8. All of the jets have 
j3 m 0.8 where we first resolve them and deceler- 
ate rapidly over a short distance in the flaring region. 
Further out, the velocities are either roughly constant 
(NGC 3 15, 0326+39 and 1553+24) or decrease more 
slowly (3C31). 
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10 20 30 

Distance from nucleus / kpc 

Figure 6: A profile of noB 1+a along the jet axis for 
the model of NGC315, deduced from the emissiv- 
ity. The energy spectrum of the radiating electrons is 
noE~ 2a+1 dE and B is the magnetic field, uq is in 
m~ 3 and B in T. Note the region of enhanced emis- 
sivity from 2.5 to 10 kpc from the nucleus. This corre- 
sponds to the knotty, quasi-helical structure in the main 
jet (Fig. 0. 



tent with flux freezing, but laminar-flow models, even 
including transverse velocity shear, do not fit quantita- 
tively (see Section^. 

As mentioned earlier, the idea that the longitudinal 
and toroidal field components are both vector-ordered 
(forming something like a simple helix) is inconsistent 
both with the observation that the brightness and po- 
larization distributions are very symmetrical under re- 
flection through their axes and with constraints from 
flux conservation. The former argument is strongest on 
large scales; weaker where the jets first brighten and 
there is significant oblique structure. The hypothesis 
that the longitudinal component has many reversals, 
but the toroidal component is well-ordered is consis- 
tent with the observations, however (our calculations 
would give identical results in this case). It is difficult 
to understand how a jet could acquire a toroidal com- 
ponent with many reversals. 

We have searched for the characteristic signature of 
Faraday rotation by a toroidal field in the jets of 3C 3 1 
and NGC 315, but see only random fluctuations which 
are most plausibly associated with foreground galaxy 
or group-scale plasma containing a disordered field. In 
any case, if our estimates of internal jet density (Sec- 
tion|4j are correct, we would expect negligible internal 



Faraday rotation even from a perfectly-ordered field. 

Although the toroidal field component is probably not 
required for the collimation of FR I jets in the regions 
we model (see below), it may be what remains of 
a dynamically-important collimating field on smaller 
scales. 

3.5 Deceleration physics 

Any theory of jet deceleration has to explain the fol- 
lowing results: 

1. Jets have (at least) two regions differentiated by 
their collimation properties: flaring and outer. 

2. The onset of rapid deceleration occurs within the 
flaring region and is quite sudden. 

3. After the rapid deceleration, three of the sources 
(0326+39, 1553+24 and NGC 315) have roughly 
constant velocities. 3C31 continues to deceler- 
ate, but more slowly. The difference may be con- 
nected to the opening angle of the conical outer 
region: 3C31 expands much more rapidly than 
the other three sources. 

4. The sudden brightening of the jets occurs within 
the flaring region and significantly before the 
rapid deceleration. It is a consequence either of 
a real increase in rest-frame emissivity or of an 
abrupt slowing of the emissivity decline. Com- 
plex and clearly non-axisymmetric structures pro- 
duce the excess emissivity. 

5. This region of enhanced emissivity persists 
roughly until the start of rapid deceleration. 
Thereafter, the jets appear somewhat smoother 
and the emissivity declines monotonically. 

6. There is evidence from the field structure of 
3C 3 1 for interaction with the external medium 



at the edge of the flaring region dLaing & Bridle, 
2002ah . but it is not clear whether this effect is 



general. 

The reaso n for the sudden brightening of the jets is un- 
clear, but Lains & Bridle ( 2002a. b) favour areconfine- 
ment shock (lSandersL ll983). as this provides a natural 
mechanism both for a sudden brightening at the same 
distance from the nucleus in both jets and for the in- 
ferred over-pressure in 3C 31 (see Section HJ). One in- 
teresting possibility is that entrainment of the external 
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medium, in the form of relatively large clumps, begins 
at the reconfinement shock, but that deceleration does 
not become significant until the entrained material has 
time to mix with the jet flow. 

3.6 Implications for unified models 

3.6.1 Total intensity 



given by Lain g & Bridj^ ( 2002ah . The dominance of 
toroidal field at large distances from the nucleus means 
that a parallel-fiel d edge is les s likely to be observed 



when 9 « 90° dLaingl I198R Although the details 



Laing & Bridlel J2002ah show the appearance of the 
3C 3 1 model over a range of angles to the line of sight, 
and demonstrate that it would be qualitatively consis- 
tent with the appearance of the extended radio struc- 
tures of BL Lac objects for small 9. 

Conversely, we find that 1553+24 has a very small 
value of 9 « 8°. This means that its side-on counter- 
parts must be very large. The length of the main jet is at 
least 60arcsec corresponding to 388arcsec (340 kpc) 
at 9 = 60°, the median angle to the line of sight. This 
is comparable in siz e with the longest jet in the B2 
sample, in NGC 3 15 dwillis et al.Lll98lh. and far larger 
than the median f^30kpc: IParma et al.ll987h . There is 
cause for suspicion unless: (a) the side-on counterparts 
of B2 1553+24 are not members of the B2 sample or 
(b) sources in that sample have linear sizes far larger 
than previously realised. 

Canvin & Laing l ( 2004 ) examined the potential selec- 
tion effects in detail and concluded that the side-on 
counterparts of 1553+24 could have escaped identifi- 
cation, either because they were missed in the original 
survey or because their angular sizes have been greatly 
underestimated. This is because there is relatively little 
lobe emission in 1553+24, so even its low-frequency 
flux is significantly affected by Doppler beaming and 
because its outer jets would be extremely faint and dif- 
ficult to detect in published images. Morganti & Parma 
(private communication) and Ledlow & Owen (2004; 
in preparation) have made more sensitive observations 
of radio galaxies in the B2 sample using the WSRT and 
the VLA in D configuration and have shown that a sig- 
nificant fraction of them have much longer radio jets 
than have previously been reported, extending many 
100's of kpc or even further. We suggest that these 
might include the missing counterparts of 1553+24. 

3.6.2 Polarization 

Our models also require that the observed polariza- 
tion structure changes as a function of angle to the 
line of sight and predictions for the 3C 31 model are 



of this effect depend on the relative amounts of ra- 
dial and longitudinal field present (these d iffer between 
sourc es), it is interesting that 3C 449 jFeretti et all 
' l999h and PKS 1333-33 dKilleen. Bicknell & Ekersl 
1986), which are extremely symmetrical and therefore 



likely to be close to the plane of the sky, show no 
parallel-field edge; neither does 0326+39 (which has 
the largest value of 9 of any of our modelled sources). 

4 Conservation-law analysis 

4.1 Assumptions and results for 3C 31 

The free models described in the previous section tell 
us the velocity and cross-sectional area of the jets. We 
can derive the external density and pressure from X- 
ray observations, and the combination allows us to 
apply cons ervation of matter, mom entum and energy. 
For 3C 31, Eaing & Bridlel |2Q02b) showed that well- 
constrained solutions exist subject to two key assump- 
tions: 

1. $ = lie, where <I> is the energy flux (with rest- 
mass energy subtracted) and II is the momentum 
flux. This must hold quite accurately if the jets 
have decelerated from bulk Lorentz factors ph 5 
on pc scales. 

2. The jets are in pressure balance with the external 
medium in their outer regions. 

For 3C 3 1 the key derived parameters at 1 kpc from the 
nucleus (the start of the modelled region) are: 

Mass flux rj 0.0005 solar masses / year. 
Energy flux « 1.1 x 10 37 W. 
Pressure sa 1.5 x 10~ 10 Pa. 
Density «2x 10" 27 kgm- 3 . 
Mach number 1.5. 

Entrainment rate sa 1.2 x 10 10 kgkpc -1 s _1 . 

Profiles of internal and external pressure, density, 
Mach number and entrainment rate are shown in Figs0 

-M 
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Figure 7: Pressure profiles alon g the axis of 3C 31 from 
a conservation-law analysis (Laing & Bridle, 2002b). 
The full line is the internal pressure, the dashed line the 
external pressure deduced from Chandra observations 
and the dotted line the synchrotron minimum pressure. 
Note that the jet is overpressured between 1 and 3.5 kpc 
from the nucleus. 
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Figure 8: Profile of internal density for 3C31, 
as deduced from the conservation-law analysis of 
Laing & Bridlel J2002hh . 
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Figure 9: Profile of internal Mach number for 3C31, 
as deduced from the conservation-law analysis of 
Lain g & Bridlel d2002hh . 
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Figure 10: Profiles of entrainment rate for 3C 31. The 
full line is deduced from the conservation-law analysis 
and the dashed line is an estimate of mass inj ection due 
to stellar mass loss (Laing & Bridle, 2002b). 
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4.2 What are the jets made of? 

Unfortunately, our analysis cannot give a definitive an- 
swer to this question, but they do imply that the jets 
are extremely light. At the innermost modelled point, 
we infer a density of 2.3 x 10~ 27 kg m~ 3 , equivalent to 
1.4 protons m~ 3 . For a power-law energy distribution 
of radiating electrons, each accompanied by a proton 
for charge balance, there must be a minimum Lorentz 
factor (the infamous 7 m i n ) in the range 20 - 50. Radi- 
ation from such low-energy electrons is at frequencies 
well below those observable. We cannot exclude any 
of the following alternatives (or, indeed, a mixture of 
them). 

1. Pure electron-positron plasma with an excess of 
particles over a power-law distribution at low en- 
ergies. 

2. Electron-positron plasma with a power-law distri- 
bution and a small amount of mixed-in thermal 
plasma. 

3. A power-law distribution of electrons, each ac- 
companied by a "cold" (i.e. not ultrarelativistic) 
proton and 7 m ; n f» 20 - 50. 

We have, however, estimated the mass injection into 
the jet from stellar mass loss within 1 kpc of the nu- 
cleus. This is commensurate with the mass flux we 
infer at that distance, and entirely consistent with the 
second alternative of electron-positron plasma and en- 
trained thermal material. We have also compared the 
mass injection rate expected from stars with our in- 
ferred entrainment rate on larger scales (Fig. [TOt, The 
two rates are comparable at 1 kpc, but an additional 
source of mass - presumably provided by boundary- 
layer entrainment - is clearly required at large dis- 
tances, where the stellar density falls off. 

5 Adiabatic models and flux-freezing 

To proceed beyond the purely empirical descriptions 
of the jets presented in Sections |2] and fJJ we must 
make further assumptions about processes that af- 
fect the components of the modelled emissivity - the 
relativistic particle energy spectrum and the strength 
of the magnetic field. The separation of the emis- 
sivity into its components is ill-determined unless 
inverse-Compton emission can be detected from the 
synchrotron-emitting electrons (not yet the case for any 




Ad tabetic, wHh same velocity end 
inilrel conditions as tree model 



Ad tabetic, wKh distributed 
particle injection 



Figure 11: Observed and model brightness distribu- 
tions for 3C 31, illustrating the effects of adiabatic en- 
ergy losses. From the top: (a) observed brightness 
distribution at a r esolution of 0.75 arcsec FWHM; (b) 
free model from Laing & Bridle ( 2002al) : (c) adia- 
batic model with the same velocity field and initial 
conditions (emissivity and field structure) as the free 
model; (d) adiabatic model with distributed particle 
injection between 2 a nd 6 kpc from the nucleus. See 



Laing & B ridle (2004) for further details. 
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Distance from nucleus / arcsec 

Figure 12: (a) Intensity profiles along the main jet axis 
for 3 C 31. The points show the observed X-ray emis- 
sion dHardcastleetal.[l2002h . the full line is the parti- 
cle injection function inferred from our adiabatic mod- 
els (arbitrarily normalised) and the dashed line is the 
observed radio emission, (b) Pressure profiles from 
Fig- El f° r comparison. The model curves in both pan- 
els start at 1 kpc from the nu cleus and do no t attem pt 
to represent the inner jet. See lLaing & B ridle (2004). 



FRI jet) and there are, as yet, no prescriptive theories 
for dissipative processes such as particle acceleration 
and field amplification or reco nnection in con ditions 
appropriate to FR I radio jets ( De Young l 2004 ). 

The energy-loss processes for the radiating particles 
can be quantified, however. It is inevitable that the par- 
ticles will suffer adiabatic losses as the jets expand and 
synchrotron and inverse Compton losses are likely to 
be negligible by comparison at radio frequencies. It is 
therefore worthwhile to compare the observations with 
models in which the radiating particles are accelerated 
before entering the region of interest and then lose en- 
ergy only by the adiabatic mechanism while the mag- 
netic field is frozen into and convected passively with 
the flow, which is assumed to be laminar. Following 
conventional usage, we refer to such models as adia- 
batic. 

Analytical expressions for magnetic-field strength and 
emissivity were gi ven for quasi -one -dimensional, non- 
relativistic jets bv lBurchl (1979) and lFanti et al.l (1982) 
and e xtended to the relativistic case by iBaum et al.l 
( 19971) . These formulae are not valid if there is veloc- 
ity shear in a direction perpendicular to any component 
of the field. We therefore developed a more general 
approach which allows the inclusion of more realistic 
geometries and velocity fields and calculates the ob- 
served brightness and polarization structure given pre- 
scribed initial conditio ns on the emissivity and field 
jLain g & Bridlel Eo04[) . We found that adiabatic mod- 
els give a fair description of the observed brightness 
and polarization distributions in the outer parts of the 
jets in 3C 31, but fail in the flaring region, where they 
predict a far steeper decline in emissivity than is ob- 
served. This is illustrated in Fig. ^2 where we show 
the brightness distribution predicted by an adiabatic 
model of 3C 3 1 , given the velocity field and initial con- 
ditions of the free model and starting at 1 kpc from the 
nucleus. The observed polarization distribution is also 
inconsistent in detail with flux-freezing in a laminar 
flow. 

We showed, however, that a modified adiabatic model 
can still be fitted to the total intensity of the flaring 
and outer regions if we add distributed injection of 
relativistic particles which then evolve adiabatically; 
the region where these particles must be injected is 
also one where there is independent evidence for re- 
cent particle accelerati on from the detection o f X-ray 
synchrotron radiation ( Ha rdcastle et all 120021) and of 
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Figure 13: Profiles o f noB 1+ct for 1553+24, from 
Canvin & Lai ng (2004). The units are as in Fig.|6] The 
full line shows the results from the free model and the 
dashed line the adiabatic model prediction, normalised 
to match at large distances from the nucleus. Note the 
good agreement over a factor of 20 in distance and 50 
in emissivity. 

a local over-pressure in the jet from dynamical argu- 
ments (Fig.0. The predicted brightness distribution is 
shown in the bottom panel of Fig. \H] and the inferred 
profile of particle injection is compared with that of the 
observed X-ray emission in Fig. El 
A simpler analysis of 1553+24 shows that the emis- 
sivity profile is fit surprisingly well by the quasi-one- 
dimensional adiabatic approximation at large distances 
from the nucleus (Fig.[T3l. We plan more sophisticated 
modelling to determine whether the slow outer regions 
of FRI jets can be adequately described by the adia- 
batic approximation. 

6 Conclusions 

We conclude that FRI jets are decelerating, relativistic 
flows, which we can now model quantitatively. Their 
three-dimensional distributions of velocity, emissivity 
and field ordering can be inferred by fitting to deep, 
well-resolved radio images in total intensity and lin- 
ear polarization. Application of conservation of mass, 
energy and momentum allows us to deduce the vari- 
ation of density, pressure and entrainment rate along 
the jets. Boundary-layer entrainment and mass input 
from stars are probably both important in slowing the 
jets. Adiabatic models and flux freezing do not work 
everywhere, but do describe the observations at large 
distances from the nucleus. Finally, we infer from the 



radio observations alone that fresh radiating particle 
must be injected where the jets are fast - this is pre- 
cisely where we detect X-ray synchrotron emission. 
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